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I. INTRODUCTION 



Probability theorists, physicists, biologists and financiers are already familiar with clas- 
sical and quantum statistical and geometric methods applied in various branches of science 
and economy [1-6]. We note that modeling of diffusion processes in nonhomogerneous me- 
dia and formulation of nonlinear thermodynamics in physics, or of dynamics of evolution of 
species in biology, requires a more extended geometrical background than that used in the 
theory of stochastic differential equations and diffusion processes on Riemann and Lorentz 
manifolds [7-10]. 

Our aim is to formulate the theory of diffusion processes on spaces with local anisotropy. 
As a model of such spaces we choose vector bundles on space-times provided with nonlinear 
and distinguished connections and metric structures [11,12]. Transferring our considerations 
on tangent bundles we shall formulate the theory of stochastic differential equations on 
generalized Lagrange spaces which contain as particular cases Lagrange and Finsler spaces 
[13-17]., 

The plan of the presentation is as follow: In Section II we give a brief summary of 
the geometry of locally anisotropic spaces. Section III is dedicated to the formulation of 
the theory of stochastic differential equations in vector bundle spaces. This section also 
concerns the basic aspects of stochastic calculus in curved spaces. In Section IV the heat 
equations in bundle spaces are analyzed. The nondegenerate diffusion on spaces with local 
anisotropy is defined in Section V. We shall generalize in Section VI the results of Section 
IV to the case of heat equations for distinguished tensor fields in vector bundles with (or 
not) boundary conditions. Section VII contains concluding remarks and a discussion of the 
obtained results. We present a brief introduction into the theory of stochastic differential 
equations and diffusion processes on Euclidean spaces in the Appendix. 
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II. NONLINEAR AND DISTINGUISHED CONNECTIONS AND METRIC 
STRUCTURES IN VECTOR BUNDLES 



As a preliminary to the discussion of stochastic calculus on locally anisotropic spaces 
we summarize some modern methods of differential geometry of vector bundles, generalized 
Lagrange and Finsler spaces. This Section serves the twofold purpose of establishing the 
geometrical background and deriving equations which will be used in the next two Sections. 
In general lines we follow conventions from [11,12]. 

Let £ = (E, 7r, F, Gr, M) be a locally trivial vector bundle, v-bundle, where F = 1Z m 
is a real m dimensional vector space, 1Z denotes the real number field , Gr = GL(m,7Z) is 
the group of automorphisms of 7Z m , base space M is a differentiate (class C°°) manifold 
of dimension n, dimM = n, ir : E — > M is a surjective map and the differentiate manifold 

E, dim E = n + m = q, is called the total space of v-bundle £ . We locally parametrize space 
£ by coordinates u a = (x\y a ), where x % and y a are, correspondingly, local coordinates on 
M and F. We shall use Greek indices for cumulative coordinates on £ , Latin indices 

k, ... = 0, 1, n — 1 for coordinates on M and a, b, c, ... = 1, 2, m for coordinates on 

F. Changings of coordinates (x k ,y a ) — > (x k ' ,y a ') on £ are written as 



A nonlinear connection, N-connection, in v-bundle £ is defined as a global decomposition 
into horizontal H£ and vertical V£ subbundles of the tangent bundle: 



x' 



.k' 



x k '{x k ),y a ' = M^'{x)y a , 




(2.1) 



T£ = H£® VS. 



(2.2) 



N-connection defines a corresponding covariant derivation in £ : 



V y A = Y l { 



OA 



+ N? (s,A)}s< 



(2.3) 



dx i 
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where s a are local linearly independent sections of E, A = A a s a and Y = Y' l Si is a vector 
field decomposed on local basis Sj on M. Differentiable functions Nf from (2.3), written as 
functions on x % and y a , i.e. as iVf (x, y) , are called the coefficients of N-connection. 
In v-bundle £ we can define a local frame (basis) adapted to a given N-connection: 

*■ = ^ = {S - = h = sr " = £ = 9 » = £>■ ( 2 - 4 ' 

The dual basis is written as 

5 a = 5u a = (5 i = 5x l = d l = dx\ 5 a = 5y a = dy a + Nf(x, y)dx l ). (2.5) 

By using adapted bases one introduces the algebra of tensorial distinguished fields (d- 
fields or d-tensors) on £,T = ®T q p s r , which is equivalent to the tensor algebra on v-bundle 
£d defined as 

n d :H£®V£^ £. 



An element t G T P J\ d-tensor field of type 



/ \ 

p r 



, can be written in local form as 



* = ttXbl--Z( x , y)hi ® - ® K ® djl ® - ® <* iq ® 

d ai ® ... ® <9 ar <8> ® ... ® 5 fes . 

In addition to d-tensors we can consider d-objects with various group and coordinate 
transforms adapted to a global splitting (2.2). For our further considerations the linear 
d-connection structure, defined as a linear connection D in £ conserving under parallelism 
the Whitney sum H£ © V£ associated to a fixed N-connection in £ will play an important 
role. Components of d-connection D are introduced as 

Torsion Tj^ and curvature Rp y $ of d-connection are defined in a standard manner: 

T(S y ,6p) = T2L6 a , 



where 



rpa _ pa pa , a 



(2.6) 



and 



where 



R(S T , Sj, Sp) — Rp yT 5 a , 



pa _ X pa _ X pa i p<P pa _ pip pa , pa <p 



(2.7) 



In formulas (2.6) and (2.7) we have used nonholonomy coefficients of the adapted frames 
(2.3) 



[8 a , Sp] = 8 a 5p - 5p5 a = wlpSp. 



(2.8) 



Another important geometric structure on £ which will be considered is the metric 
structure on G, defined as a nondegenerate, second order, with constant signature (in this 
work positive), tensor field G a p on £. We can associate a map 

G(u) : T u £ xT u £ 



to a metric 



G = G aP du a <g> dvP 



parametrized by a symmetric, rankG = q, matrix 



G a j G a i) J 



where = G (<%, dj) ,G ia = G (<%, d a ) ,G ab = G (d a , d b ) . 

We shall be interested by such a concordance of N-connection and metric structures when 
G (Si, d a ) = 0, or, equivalently, 



N«(x,y) = G lb (x,y)G ba (x,y), 



(2.9) 
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where matrix G ba is inverse to matrix G a b- In this case metric G on E is defined by two 



independent d-tensors, of type 
respect to bases (2.4) and (2.5) as 



\ 00 / 



and h a b of type 



\ 02 / 



, and can be written with 



G = gij(x, y)dx l <g> dx 3 + h ab (x, y)5y a ® Sy b 



We shall call a metric G compatible with a d-connection D if conditions 



(2.10) 



D a G 



Pi 







(2.11) 



are satisfied. 

Definition 1. A vector bundle £ on base M provided with compatible N-connection, 
d-connection and metric structures (when conditions (2.9) and (2.11) are satisfied) is called 
a space with local anisotropy, (la-space) and denoted as £^. 

Remarks. 

1. For the case when instead of v-bundle £ the tangent bundle TM is considered 
the conditions (2.9), (2.11) and the requirement of compatibility of compatibility of N- 
connections with the almost Hermitian structure on TM lead to the equality [11,12] (see 
metric (2.10) on TM ) 



hij(x,y) = gij(x,y). 

The metric field gij(x,y) is the most general form of metric structure with local anisotropy, 
considered in generalized Lagrange geometry. 

2. Metrics on a Lagrange space (M, C) can be considered as a particular case of metrics 
of type (2.10) on TM for which there is a differentiable function C :TM — > 1Z with the 
property that d-tensor 

1 d 2 C 



9%j{x,y) 



2 dtfdyi 



(2.12) 



is nondegenerate. 



3. We obtain a model of Finsler space (M, L) if C = L 2 , where L is a Finsler metric 
[11-12]. 

We emphasize that in this work all geometric constructions and results on stochastic 
calculus will be formulated for the general case of la-spaces. 

On a la-space E^ we can consider arbitrary compatible with metric G d-connections 
Fg , which are analogous of the affine connections on locally isotropic spaces (with or not 
torsion). On E^ it is defined the canonical d-structure T ^ with coefficients generated by 
components of metric and N-connection 



jk = Ljk-> r j a = Cj a , r a j = o, r ab = o, 

~?]k = 0, ~?] b = 0, K = L a bk , K = (2.13) 

where 

1 . 

L 'jk = ^9 W {hg P i + bj9 P k - $p9jk), 

L a bi = d b N? + h^hc - d b N?h dc - d c N?h db ), 

C) c = \g ik d b9jk , 

C a bc = \h ad (d c h db + d b h dc - d d h bc ). 

In formulas (2.13) we have used matrices and h ab which are respectively inverse to matrices 
gij and h ab . 

We also present the explicit formulas for unholonomy coefficients u^ 7 of the adapted 
frame basis (2.4): 

™ k H = 0, w k aj = 0, w k ia = 0, w k ab = 0, w c ab = 0, 

K = R h wb ai = ~ d « N l < = 9 « N l (2-14) 

Putting (2.13) and (2.14) into, correspondingly, (2.6) and (2.7) we can computer the com- 
ponents of canonical torsion T ^ and curvature R ^ s with respect to the locally adapted 
bases (2.4) and (2.5) (see [10,11]). 
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Really, on every la-space £V a linear multiconnection d-structure is defined. We can 
consider at the same time some ways of local transports of d-tensors by using, for instance, 
an arbitrary d-connection Pg , the canonical one r ^ , or the so-called Christoffel d-symbols 
defined as 

{^} = l -G aT {5,Gp T + SpGyr - 5 T G Pl ). (2.15) 

Every compatible with metric d-connection can be characterized by a corresponding 
deformation d-tensor with respect, for simplicity, to {j^} '■ 

= (2-16) 

(the deformation of the canonical d-connection is written as 

(2.17) 

Perhaps, it is more convenient to consider physical models and geometric constructions with 
respect to the torsionless d-connection {-^}.The more general ones will be obtained by 
using deformations of connections of type (2.16). But sometimes it is possible to write out 
d-covariant equations on £^ , having significance, by changing respectively components {^} 
on . This holds for definition of stochastic differential equations on la-spaces (see, in 
particular, [1,33] on diffusion on Finsler spaces) and in this paper we shall use the last way. 

Let suppose that £j\f is locally trivial and a-compact. In this case £jy is a paracom- 
pact manifold and has a countable open base. We denote as F{£^- ) the set of all real 
C°°-functions on and as F (£jv) the subclass of F (£j\f) consisting from functions with 
compact carriers. F (£^) and F (£^f) are algebras on the field of real numbers 1Z with 
usual operations / + q, fq and A/(/, q G F(£j^ ) or F (£j^) ,\&1Z). 

Vector fields on £^ are defined as maps 

V : u e £ M ^V(u)eT u (£ M ) . 

Vectors (d a ) u , (a = 0, 1, 2, m + n — 1), form a local linear basis in T u (£^) . We shall also 
use decompositions on locally adapted basis (2.4), (5 a ) , and denote by X(^) the set of 
C°°-vector fields on £^. 

8 



Now, let introduce the bundle of linear adapted frames GL (£m) on £-N ■ As a linear 
adapted frame e = [e^], (a = 0, 1, m+n— 1), in point u G £m we mean a linear independent 
system of vectors e„ G T u (£V) obtained as a linear distinguished transform of local adapted 
basis (2.4), i.e.e a = e£5 a , where G GL d (72) = GL (n, 72 ) ® GL (m, 72) . Then GL (£ N ) is 
defined as the set of all linear adapted frames in all points u G £/v : 

GL (£/v0 = = (m, e), m G £\f an d e is a linear adapted frame in the point u}. 

Local coordinate carts on GL(£x) are defined as (u a ,<Pa) , where U a = {r = 
(u, e) G GL (£j\f) , m G U a C £V, and e is a linear adapted frame in the point u} , 
^ (r) = (ip a (u) = (u a ), e^j and e« = e^<5 Q \ u . So GL (£/y) has the structure of G°°-manifold 
of dimension n + m + n 2 + m 2 . Elements a G GL rf (72.) act on GL (£/v) according the for- 
mula T a (w, e) = (u, ea), where (ea)„ = aie^. The surjective projection 7r : GL (£/v) — > £a/" is 
defined in a usual manner by the equality 7t(m, e) = w. 

Every vector field L G X (£/v) induces a vector field L on GL {Em) ■ Really, for / G X (£/v) 
we can consider 

(If) (r) = ^/((exptL)u,(exptL),e) | t=0 , (2.18) 

where r = (u, e) and 

(exptL)*e = [(expiL)*ei, (exp£L)*e 2 , (exp£L)*e m+n ], 

is the differential (an isomorphism T u (E^- ) — > Ti ey . ptL \E^f f° r every u G £V) of exptL and 
the local diffeomorphism u — > v(t,u) is defined by differential equations 

— (t,u) = a a (v(t,u)), (2.19) 

(L = a a {u)5 a ) , v (0, w) = -u. 

Let L G X and introduce functions f"(r) G L (GL(£^)) for every a = 0, 1, 2, m+ 
n — 1 by the equalities 

= (e- 1 ) V(u) (2.20) 



written in locally adapted coordinates r = (u a , e = (e™)) on the manifold GL (£tf) , where 
L = a a {u)5 a and e~ l is the matrix inverse to e. Because the equality (2.20)does not depend 
on local coordinates, we have defined a global function on GL (£j\r) ■ It's obvious that for 
L(i), L(2) G X w e have 

{ki)fL m ) (r) = /[L (1) L (2) ] W, 
where L(i) and Z( 2 ) are constructed similarly to operator (2.18), and 

[L(i),L( 2 )] = ^(1)^(2) - L(2)L {1) . 

A distinguished connection defines the covariant derivation of d-tensors in £^ in a usual 
manner. For example, we can introduce a d-covariant derivation DB of a d-tensor field 
5 («) = («) in the form 



E r^j («) b%£^{u) - f („) 3,X:&>), (2-21) 

e=l r=l 

or the covariant derivative LV-B in the direction F = G X (£yy ) , 

and the parallel transport along a (piecewise) smooth curve c : 1Z Dl = (tit 2 ) 3 t — > c(t) 
(considering B (t) — B (c(t)) ) 



j t BZ:.t(t) + E r«! (c(0) «;i;^(cW)^- 



E ^ (<#)) ^X:& 9 (c(0)^ = 0. (2.23) 
For every r G GL (£c) we can define the horizontal subspace 

H r = {U = a a 5 a \ u -Y a M (u) e>^, a a G TZ m+n } 

of T u (GL(£x)). Vector U G H r is called horizontal. Let £ G T u (£/v ) , then £ G 
T r (GL(£tf)) is a horizontal lift of £ if the vector £ is horizontal, i.e. ir (r) = u and 
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(dir) r £ — £. If r is given as to satisfy ir (r) = -u the £~ is uniquely defined. So, for given 
[/ G X(£_a/) there is a unique [/ G X (GL (£V)) , where U r is the horizontal lift U n (r) for 
every r G GL (£jy) XJ is called the horizontal lift of vector field U. In local coordinates 
U = U a (u) 5 a - F 5 ap (u) U a (u) ef ^ if U = U a (u) S a . 
In a similar manner we can define the horizontal lift 

c(t) = (c(t),e(t)) = [e (t),e 1 (t),...,e^ 1 (t)} G GL (£ M ) 

of a curve c (t) G Em with the property that n (c(t)) = c(t) for t G / and ^(t) is horizontal. 
For every a = 0,1, ...q — 1 there is a unique vector field L G X (GL(Sx)) , for which (^«) is 
the horizontal lift of vector G T u (^V) for every r — (u, e — [e , e\, e g _i]) . In coordinates 
we can express 

z a = - a4 e ?4r (2 - 24) 

Vector fields form the system of canonical horizontal vector fields. 

Let B (u) = B^p*'''p* (u) be a (p,s)-tensor field and define a system of smooth functions 

•B l/J — l r J — ^1(52...^ l M J e 7i e 72---e7 P e^e^...e^| 

(the scalarization of the d-tensor field .£> (r) with the respect to the locally adapted basis e) 
on GL ^ where we consider that 

0102-«p / n 01 02 /3s 



, \ AXJ^U./ .. .Lip , . 



• ••e~ <x) e— <x> e— ... <K> e 



* 5 



the matrix is inverse to the matrix e^, basis e* is dual to e and r = (w, e). It is easy to 
verify that 

^411) (0 = (^VB)^f...^ (r) (2.25) 
where the covariant derivation V^A- = A§i is taken by using connection 

r- — pOLpPpir 4- p 7 p-A p°" 

1 7 — e a e e 7 /3 7 ^ e 7 e o-°7 e /3 
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in GL (£j\r) (induced from Sj^f). 

In our further considerations we shall also use the bundle of orthonormalized adapted 
frames on £jy, defined as a subbundle of GL (Sj^) satisfying conditions: 

O (£n) = {r = (u, e) G GL (£jy) , e is a orthonormalized basis in T u (£m)}- 

So r = (u a , G O {£ N )) if and only if 

Gapele^ = (2.26) 

or, equivalently, 

where the matrix G a/3 is inverse to the matrix G a p from (2.10). 

III. STOCHASTIC DIFFERENTIAL EQUATIONS IN VECTOR BUNDLE 

SPACES 

In this Section we assume that the reader is familiar with the concepts and basic results 
on stochastic calculus, Brownian motion and diffusion processes (an excellent presentation 
can be found in [7-9, 18-20], see also a brief introduction into the mentioned subjects in the 
Appendix of this paper). The purpose of the Section is to extend the theory of stochastic dif- 
ferential equations on Riemannian spaces [7-9] to the case of spaces with general anisotropy, 
defined in the previous Section as v-bundles. 

Let A~, A~, A~ G X (£m) an d consider stochastic differential equations 

dU(t) = A~odB«(t) + A d (U(t))dt, (3.1) 

where a = 1, 2, r and o is the symmetric Q-product (see Appendix Al). We shall use the 
point compactification of space £^ and write £ n = £ or £ ^ = £ U{A} in dependence 
of that if £j^f is compact or noncompact. By W (£tf) we denote the space of paths in £^ , 
defined as 

W {£m) — { w '■ w is a smooth map [0,oo) — > £ n with the property that w (0) G £jv and 
w(t) = A, w(f) = A for all t' > t} 
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and by B [W (E M )) the cr-field generated by Borel cylindrical sets. 
The explosion moment e(w) is defined as 

e(w) = M{t,w(t) = A} 

Definition 3.1. The solution U = U (t) of equation (3.1) in v-bundle space E^ is defined 
as such a {Ft) -compatible W (£/v)- valued random element (i.e. as a smooth process in Ej\f 
with the trap A), given on the probability space with filtration (Ft) and r- dimensional 
Tt)-Brownian motion B = B(t), with B(0) = 0, for which 

f{U(t))-f{U(0)) = 
f A~ (t) (U (s)) 5B« (s) + f (A o f) (U (s)) ds (3.2) 

JO JO 

for every f G Fo(£j\f) (we consider / (A) = 0), where the first term is understood as a 
Fisk-Stratonovich integral. 

In (3.2) we use 5 B a (s) instead of dB a (s) because on the Brownian motion must be 
adapted to the N-connection structure. 

In a manner similar to that for stochastic equations on Riemannian spaces [7] we can 
construct the unique strong solution to the equations (3.1). To do this we have to use the 
space of paths in W starting in point 0, denoted as Wq, the Wiener measure P w on Wq, cr- 
field B t (VF^-generated by Borel cylindrical sets up to moment t and the similarly defined 
cr-field. 

Theorem 3.1. There is a function F : x -> W (E/j- ) being r\^B(£ N ) x B t 
(W£T xpW /B t (W (Ft)) -measurable (index ji runs all probabilities in (E^, B(Ej^) ) ) for 
every t > and having properties: 

1) For every U(t) and Brownian motion B = B (t) the equality U = F{U (0) ,B) a.s. is 
satisfied. 

2) For every r- dimensional [Ft)-Brownian motion B = B (t) with B — B (0) , defined 
on the probability space with filtration Ft, and E^-valued F,-measurable random element £, 
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the function U = F(£,B) is the solution of the differential equation (3.1) with U (0) = £, 
a.s. 

Sketch of the proof. Let take a compact coordinate vicinity Vwith respect to a locally 
adapted basis 5 a and express Ag_ = a" {u) 5 ai where functions cr" {u) are considered as 
bounded smooth functions in Jl m + n ) and consider on V the stochastic differential equation 

dU? = a* (Un o SB" (t) + < (U t ) dt, (3.3) 

[tf = «",(« = 0,1, 
Equations (3.3) are equivalent to 

dU? = a%(U?)dB°(t) + a%(U t )dt, 

1 



where a% (u) = a% (u) + \ ££ =1 [~i^J («) • It's known [7] that (3.3) has a unique 
strong solution F : K n+m x W -> W n+m or F(«, w) = (U(t, u,w)) . Taking t v (w) = m£{t : 
C/(t, G V } we define 

C/v (t, u,w) = U(tf\ r v (u>) , «, w) (3.4) 

In a point u E V fl V, where V is covered by local coordinates u a , we have to consider 
t r ansf or mat ions 

^(5(«))=^(«)|£, 

where coordinate transforms w a (-u a ) satisfy the properties (2.1). The global solution of (3.3) 
can be constructed by gluing together functions (3.4) defined on corresponding coordinate 
regions covering E^- . (} (end of the proof) 

Let P u be a probability law on W of a solution U — U (t) of equation (3.1) with 
initial conditions C/(o) = u. Taking into account the uniqueness of the mentioned solution 
we can prove that U — U (t) is a A-diffusion and satisfy the Markov property [7] (see also 
Appendix A3). Really, because for every / e F {Em) 
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df (U (t)) = (A~f) (U (t)) o dw" + (A f) (U (t)) dt = 
(A~f) (U (t)) du? + (A f) (U (t))d+ \d (A~f) (U (t)) ■ duf it) 

and 

d (Aff) (U (t)) = A Q (Atf) (U (t)) o duf (t) + (A^f) (U (t)) dt, 

we have 

d (A-f) (U (t)) ■ dw« (t) = ±A^ (A~f) (U (t)) dt. 
Consequently, it follows that 

df (U (t)) = (A~f) (U (t)) du? (t) + (Af) (U (t)) dt, 

i.e. the operator (Af) , defined by the equality 

Af= 1 -t^(A Q f) + A f (3.5) 

a=l 

generates a diffusion process {P u },u G Ej\f. 

The above presented results are summarized in this form: 

Theorem 3.2. A second order differential operator Af generates a A-diffusion on 
W(£j^) of a solution U = U (t) of the equation (3.5) with initial condition U (0) = u. 

Using similar considerations as in flat spaces [7] on carts covering E^,, we can prove the 
uniqueness of A- diffusion {P„}, u G £/v on W(Ej^). 

IV. HEAT EQUATIONS IN BUNDLE SPACES AND FLOWS OF 

DIFFEOMORFISMS 

Let v-bundle E^ be a compact manifold of class C°°. We consider operators 

A ,A 1 ,...,A r G X (Ex) 
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and suppose that the property 

E[Supte[o,i]Sup ueU \D-{f(U (t,u,w))}\] < oo 

is satisfied for all f E F (£n) and every multiindex a in the coordinate vicinity U with U 
being compact for every T > 0. The heat equation in F (£n) at is written as 

dv 

—{t,u)=Av{t,u), (4.1) 

lim v(t,u) = f(u), 

where operator A acting on F(F (£jv)) is defined in (3.5). 

We denote by C 1,2 ([0, oo) x F (£j^)) the set of all functions f(t, u) on [0, oo) x £^ being 
smoothly differentiable on t and twice differentiable on u. 

The existence and properties of solutions of equations (4.1) are stated according the 
theorem: 

Theorem 4.1. The function 

C(t, u) = E[f(U(t, u, w)} G C°° QO.oo) x £ N ) (4.2) 

/ G F (S^f) satisfies heat equation (4-1)- Inversely, if a bounded function u{t,u) G 
C 1 ' 2 ([0,oo) x £tf) solves equation (4-1) and satisfies the condition 

lim E[u(t — dfc, U (ufc, u, w)) o~k < t] = (4.3) 

for every t > and u G £j^, where a k = inf{t, U(t, u, w) G D k } and D k is an increasing 

sequence with respect to closed sets in £^, \JD k = £^. 

k 

Sketch of the proof. The function ((t,u) is a function on £jy, because u — > 
/ (U (t, u, w)) G C°° and in this case the derivation under mathematical expectation symbol 
is possible. According to (3.2) we have 

/ (U (t, u, w)) — f (u) — martingale + / (At) (U (s, u, w)) ds 

Jo 

for every u G £_m, i.e. 
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C(t, u) = f (u) + f E[(Af) (U (s, u, w)) ds. (4.4) 



Because A n f G F (8^) , (n — 1, 2, ...) , we can write 

C(t, u) = f(u)+t (Af) (u) + J* dt x jf* E [(A 2 f) (U (t 2 , u, w))] dt 2 = 

f(u) + t (Af) (u) + ^ (A 2 f) (u) + jT dix jT* 1 rft 2 ^ £[(A 3 /) (C/ (*s, «, w))]df 3 = 



/ («) + * (Af) (U) + ^ (A 2 f) (U) + ... + jf* d*x jT* 1 d* 2 ... jf'"" (C/ «, ™)) < 

from which it is clear that 

C(t,u)eC°°([0,oo)xS M ). 
In Chapter V, Section 3 of the monograph [7] it is proved the equality 

(A( t )(u)=E[(Af)U(t,u,w))} (4.5) 

for every t > 0, where ( M ) = C (t, u) ■ 
From (4.4) and (4.5) one follows that 



C(t,u) = f(u)+ f(AQ (s,u)ds 
Jo 



and 

^(t,u) = A( (t,u), 

i.e. C = C (t, u ) satisfies the heat equation (4,1). 

Inversely, let v(t,u) G C 1 ' 2 ([0, oo) x £j^) be a bounded solution of the equation (4.1). 
Taking into account that P (e[U (-,u,w))] = oo) = 1 for every u G £n and using the Ito 
formula (see (A6) in the Appendix) we obtain that for every t > and < t < t . 

E\v (t -t f\a n ,U (t f\a n ,u,w^} -v(t ,u) = 

ft y\ <T„ Qjy 

E[J o {(Au)(t -s,U(s,u,w))- — (t -s,U(s,u,wj)}ds}. (4.6) 
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Supposing that conditions (4.3) are satisfied and considering nfoowe obtain 
E = {v (t — t,U (t, u, w)) ; e[U (•, u, w)] > t} = v (t , u) . 

For t t we have E[f (U (t , u, w))] = ( (t , u) , i.e. v (t, u) = ( (t, u) .<> 

Remarks; 1. The conditions (4-3) are necessary in order in order to select a unique 

solution of (4-1) ■ 
2. Defining 

C (t, u) = E[exp{ f C (U (s, u, w)) ds}f (U (t, u, w))} 
Jo 

instead of (4-2) we generate the solution of the generalized heat equation in £j\f : 

dv 

— (t, u) = (Au) (t, u) + C (u) v (t, u) , 
lim v (t,u) = f (u) . 

t[Q,u— >u 

For given vector fields e X(^ A r), (a) = 0, 1, ...,r in Section III we have constructed 
the map 

U = (U(t,u,w)) :S M xW r 3 (u,w) -+U{-,u,w) G W{£ N ) , 
which can be constructed as a map of type 

[0, oo) x £m x Wq 3 (u, w) — > U (t, u, w) e£u . 

Let us show that map u G £j\f — > U (t, u, w) G £m is a local diffeomorphism of the manifold 
£tf for every fixed t > and almost every w that G £j^ . 

We first consider the case when £ M = TZ n+m , a (u) = (o% («)) G K n+m <g> U n+m and 
b (u) = (b a (u)) G TZ n+m are given smooth functions (i.e. C°°-functions) on TZ n+m ) || a (u) \\ 
+ || b (u) ||< K (1 + |w|) for a constant > and all derivations of o a and 6 a are bounded. 
It is known [7] that there is a unique solution U = U (t, u, w) , with the property that 
E[(U (t)) p ] < oo for all p > 1, of the equation 

dU? = (U t ) du? (t) + b a (U t ) dt, (4.7) 

Uq = u,(a = 1,2, ...,m + n — 1), 
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defined on the space (w o r , P w ) with the flow (JF°) . 

In order to show that the map u — > U (t, u, w) is a diffeomorphism of TZ n+m it is more 
convenient to use the Fisk-Stratonovich differential (see, for example, the Appendix Al) and 
to write the equation (4.7) equivalently as 

dU t a = al (U t ) o 5uf it) + b a (U t ) dt, (4.8) 

U = u, 



by considering that 



r («) = & a ( M ) + Jt °£ («) • (4.9) 



2.. 

a=l 



We emphasize that for solutions of equations of type (4.8) one holds the usual derivation 
rules as in mathematical analysis. 

Let introduce matrices o'~ = (a' (u)|g = ^crS ( M )) , b' (u) = (if (u) a p = g^) ,1 = 8% and 
the Jacobi matrix Y (t) = (Y^ (t) = j^- (t, u, wfj , which satisfy the matrix equation 

Y(t)=I+ f a' Tv (U (s)) Y (s) o du? (s) + /' b' (U (s)) Y (s) ds. (4.10) 

JO JO 

As a modification of a process U (t, u, w) one means a such process U (t, u, w) that 
P W {U (t, u,w) = U (t, u, w) for all t > 0} = 1 a.s. 

It is known this result for flows of diffeomorphisms of flat spaces [24-26,7]: 

Theorem 4.2. Let U (t,u,w) be the solution of the equation (4-8) (or (4-7)) on Wiener 
space (Wq, P w ^j . Then we can choose a modification U (t, u, w)of this solution when the map 
u — > U (t,u,w) is a diffeomorphism TZ n+m a.s. for every t e [0, oo). 

Process u — U (t, u, w) is constructed by using equations 

dU? = aZ(U t )o5v?(t)-b a (U t )dt, 
U = u. 

Then for every fixed T > we have 

U(T -t,u,w) = U (t, U (T, u, w) , w) 
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for every 0< t < T and u P w -a.s., where the Wiener process w is defined as w (t) = 
w (T - t) - w (T) ,0 < t < T. 

Now we can extend the results on flows of diffeomorphisms of stochastic processes to 
v-bundles. The solution U (t,u,w) of the equation (3.1) can be considered as the set of 
maps U t : u — > U (t, u, w) from £^ to ~E N = £^ U {A}. 

Theorem 4.3. A process \U\ (t,u,w) has such a modification, for simplicity let denote 
it also as U (t, u, w) , that the map U t (w) : u — > U (t,u, w) belongs to the class C°° for every 
f G F (£jv) and all fixed t G [0, oo) a.s. In addition, for every u G U and t G [0, oo) the 
differential of map u — > U (t, u, w) , 

U (t,u,w), : T u (U (t,u,w)) -> T u{tjU:W) (£ M ) , 

is an isomorphism, a.s., in the set {w : U (t,u,w) G £j^ }. 

Proof. Let u G £m and fix t G [0, oo) . We can find a sequence of coordinate carts 
Ui, Uii ---jUp C £m that for almost all w that U (t, u C w) G £a/ there is an integer p > 
that {U (s,u ,w) : s G [(k — l)t/p,kt/p]} C U\\ ,(k = 1,2, ...p). According to the theorem 
4.2 we can conclude that for every coordinate cart U and {U (s,u ,w) ; s G [0, 1]} C U a 
map v — > U (to, v, w) is a diffeomorphism in the neighborhood of t>o.The proof of the theorem 
follows from the relation U (t, wo, w) = [U t / P (8( p -i)t/ p u?J o ... o U t / P (6t/pW^j o U t / P ] (u ) , where 
6 t : W'q -> W r Q is defined as (9 t w) (s)=w(t + s)-w (t) .<> 

Let Aq, Ai, A r G X (£j\r) and U t = (U (t,u,w)) is a flow of diffeomorphisms 
on £/v . Then A§, A\, A r G ~X(GL (£//)) define a flow of diffeomorphisms r t = 
(r(t,r, w)) on GL(£j^) with (r (i, r, it;)) = ([/ (t, w, it;) , e (t, w, it;)) , where r = (u, e) and 
e(t,r,u) = U (t, u, w)^ e is the differential of the map u — > U (t,u,w) satisfying the 
property U(t,u,w)^e = [U (t, u, w)^ eo, U (t, u, w)^ e±, U (t, u, w)^ e ? _i]. In local coordi- 
nates A~ (u) = o-~5 a , (a = 1, 2, r) , A (u) = b a (u) 5 a , (t, u, w) = (t, u, w) e^, where 

(t, u,w) is defined from (4.10). So we can construct flows of diffeomorphisms of the 
bundle £^ . 
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V. NONDEGENERATE DIFFUSION IN BUNDLE SPACES 



Let a v-bundle E^ be provided with a positively defined metric of type (2.10) being 

compatible with a d- connection D = {Fg 7 }. The connection D allows us to roll £/v along 

a curve 7 (t) C 'R n+m in order to draw the curve c (t) on E^ as the trace of 7 (t) . More 

exactly, let 7 : [0, 00) 3 t -> 7 (t) C ft n+m be a smooth curve in TZ n+m , r = (u,e)eO (Ex) ■ 

We define a curve c (£) = (c (£) , e (£)) in O (£//) by using the equalities 

dc a (t) _ drf± 
dt ~ 6 ^ [ ) dt ' 

(5.1) 

c 1 >(0)=„",ej(0)=ej. 

Equations (5.1) can be written as 

d ^ = l a (c(t))di\ 
c(0)=r, 

where{I/ a } is the system of canonical horizontal vector fields (see (2.21)). Curve c(t) = 
n (c (£)) on £jv depends on fixing of the initial frame p in a point u; this curve is parametrized 
as c(t) = c (t, r, 7) , r = r (u, e) . 

Let u> (£) = (w-(t)) is the canonical realization of a n+m-dimensional Wiener process. 
We can define the random curve U (t) C Em in a similar manner. Consider r (£) = (r (£, r, u>)) 
as the solution of stochastic differential equations 

dr(t) =La(r(t))o5w*(t), (5.2) 

r (0) = r, 

where r (t, r, iu) is the flow of diffeomorphisms on O (Ejy ) corresponding to the canonical 
horizontal vector fields Li,L 2 , and vanishing drift field L = 0. In local coordinates 

the equations (5.2) are written as 

dU a (t)=el(t)o5w^(t), 
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deZ(t) = -r% r (U(t))elo6uf>, 

where r (t) = (u a (t) , (f)) . It is obvious that r (t) = ([/ Q (f ) ,e«(t))eO (£^ ) if r (0) G 
O (<£V ) because Lg_ are vector fields on O (Ej^ ) . The random curve {U a (t)} on Ex is 
defined as U (t) = n [r (t)] . We point out that aw = (aw (t)) is another (n+m)-dimensional 
Wiener process and as a consequence the probability law U (-,r,w) does not depend on 
a E O (n + m) . It depends only on u = n (r) . This law is denoted as P w and should be 
mentioned that it is a Markov process because a similar property has r (•, r, w) . 

Remark 4.1. We can define r (t,r,w) as a flow of diffeomorphisms on GL (Ej^ ) for 
every d-connection on Em- In this case 7r [r (•, r, w)] does not depend only onu = n (t) and in 
consequence we do not obtain a Markov process by projecting on E^. The Markov property 
of diffusion processes on Ejy is assumed by the conditions of compatibility of metric and 
d-connection (2.11) and of vanishing of torsion. 

Now let us show that a diffusion {P u } on E^ can be considered as an A-diffusion process 
with the differential operator 

A=^A e +b, (5.3) 
where is the Laplace-Beltrami operator on E^f, 

A £ /^ZW^^-{^, (5.4) 

where operator D a is constructed by using Christoffel d-symbols (2.15) and b is the vector 
d-field with components 

Theorem 5.1. The solution of stochastic differential equation (5.2) on O ) defines 
a flow of diffeomorphisms r (t) = (r (t, r, w)) on O (Ej^ ) and its projection U (t) = n (r (t)) 
defines a diffusion process on Ej^ corresponding to the differential operator (5.3). 

Proof. Considering / (r) = / (u) for r = (u, e) we obtain 

/ (U(t)) -f(U (0)) = / (r (0) - / (r (0)) = jT (Z,/) (r ( a )) o 8vP = 
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f L a f (r (a)) ^ + l f g La ( r ( s) ) dg. 

Let us show that = X) L a (L a f) = Af. Really, because the operator (5.3) can be written as 

and taking into account (2.25) we have 

4 (4/) = 4 (*W) a = (Fw/U = (V 7 V,/) elel 

Now we can write 

E 4 (4/) = A/ )«£ef = G^D a D p f 

(see (2.26)), which complete our proof. ^> 

Definition 5.1. The process r (t) = (r(t,r,w)) from the theorem 5.1 is called the 
horizontal lift of the A-diffusion U (t) on E^. 

Proposition 5.1. For every d-vector field b = b a (u) 5 a on E^ provided with the canon- 
ical d-connection structure there is a d-connection D = {T^ 7 } on E^ , compatible with 
d-metric G a/ 3,which satisfies the equality (5.5). 

Proof. Let define 

where b a = G a pb^ . By straightforward calculations we can verify that d-connection (5.6) 
satisfies the metricity conditions 

and that 

We note that a similar proposition is proved in [7] for, respectively, metric and affine 
connections on Riemannian and affine connected manifolds: M. Anastasiei proposed [28] 
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to define Laplace-Beltrami operator (5.4) by using the canonical d-connection (2.13) in 
generalized Lagrange spaces. Taking into account (2.16) and (2.17) and a corresponding 
redefinition of components of d- vector fields (5.5), because of the existence of multiconnection 
structure on the space H 2n , we conclude that we can equivalently formulate the theory of 
d-diffusion on if 2ri -space by using both variants of Christoffel d-symbols and canonical d- 
connection. 

Definition 5.2. For A = |Af an A-diffusion U (t) is called a Riemannian motion on 
Let an A-differential operator on Ex is expressed locally as 

where / e F (£jv) > matrix a af3 is symmetric and nonegatively defined . If a a/3 (u) > for 
all u and £ = (£ a ) e lZ q \{0}, than the operator A is nondegenerate and the corresponding 
diffusion is called nondegenerate. 

By using a vector d-field b a we can define the 1-form 

^(6) = b a (u) 5u a , 

where b = b a 5 a and b a = G^b 13 in local coordinates. According the de Rham-Codaira 
theorem [27] we can write 

uq,) = dF + 5(3 + a (5.7) 

where F G F (£/v) , f3 is a 2-form and a is a harmonic 1-form. The scalar product of p-forms 
(^V ) on is introduced as 

(o£,(3)b= / <a,(3>5u, 

where 

a = ^ a 7l72 ... 7p <5w 71 /\ <5m 72 A ••• A 5m> > 

7i<72<---<7p 

/? = E /3 7172 ... 7p ^ 7l A^ 72 A-A^ 7p > 

71<72<---<7p 
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^7l72---7p _ Q11T1Q12T2 Q^P T pp 



TIT2---Tp ) 



<a,(3>= J2 «7 1 7 2 ...7 P («)/3 7172 - > 

7i<72<---<7p 



8u = yj\ det G a(3 \Su°Su 1 ...Su q - 1 . 
The operator 8 : A p (£^) ~^ \-i {^n) from (5.7) is defined by the equality 

(da, P) p = (a, 6P) , a e A p _i {£ N ) ,/3eA p {£ N ) ■ 

De Rham-Codaira Laplacian □ : A p (£x) ~^ h-p {^n) is defined by the equality 

U = -(S+5d). (5.8) 

A form a e A p (£/y) is called as harmonic if Da = . It is known that Da = if and 
only if da — and 5a = O.For / e F (£j\r) and U G X we can define the operators 
gradf G X (£V) an d <iit>£/ G F (£V) by using correspondingly the equalities 

<?™# = W (5.9) 

and 

diuC/ = -Su v = . 1 £ a f^yidetGl) . (5.10) 
^/| det C) V / 

The Laplace-Beltrami operator (5.3) can be also written as 

A e f = div(gradf) = -66 f (5.11) 

for F (M) . 

Let suggest that £j\f is compact and oriented and {P u } be the system of diffusion measures 
defined by a A-operator (5.3). Because £^ is compact P u is the probability measure on the 
set W (£tf) = W {£*[) of all continuous paths in £^ . 

Definition 5.3. The transition semigroup T t of A- diffusion is defined by the equality 

(T t f)(u)= J f{w{t))P u {dw),feC{£ N ). 
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For a connected open region Q C £/v we define p n w G W (Q) , w G W (Ex) by the 
equality 

IP {V ~ \A,if t> m (w), 

where tq (w) = inf{t : w (t) ^ Q}. We denote the image-measure P u (w G Q) on map p n as 
P^ ; this way we define a probability measure on W which will be called as the minimal 
A-diffusion on Q. The transition group of this diffusion is introduced as 



(??/)(«) = / = 
J f(w (t)) I {m(w)>t} P u (dw) JeC p (tt). 



Definition 5.4. The Borel measure \x (du) on E^ is called an invariant measure on 
A-diffusion {P u } if J T t f (u) /i (du) = J f (u) fj, (du) for all f G C (£ A f) . 

Definition 5.5. An A-diffusion {P u } is called symmetrizable (locally symmetrizable) if 
there is a Borel measure v (du) on E^ {v n (du) on Vtj that 

J T t f(u)g(u)u(du) = J f(u)T t g(u)u(du) (5.12) 
for all f,g G C (Em) an d 

( J Tff (u) g (u) u n (du) =//(«) T?g (u) v Q (du) 
n n 

for all f,geC(Q)). 

The fundamental properties of A-diffusion measures are satisfied by the following theorem 
and corollary: 

Theorem 5.2. a) An A-diffusion is symmetrizable if and only if 5(3 = a = (see (5.8)); 
this condition is equivalent to the condition that b = gradF, F G F (Ej^) and in this case the 
invariant measures are of type Cexp[2F (u)]du, where C = const. 

b) An A-diffusion is locally symmetrizable if and only if 5(3 = (see (5.8)) or, equivalently 
, dw(p) = 0. 
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c) A measure cdu (constant c > 0) is an invariant measure of an A-diffusion if and only 
if dF = (see (5.8)) or, equivalently, 5wi b \ = —divb = 0. 

Corollary 5.1. An A-diffusion is symmetric with respect to a Riemannian volume du 
(i.e. is symmetrizable and the measure v in (5.14) coincides with du) if and only if it is a 
Brownian motion on E^. 

We omit the proofs of the theorem 5.2 and corollary 5.1 because they are similar to those 
presented in [7] for Riemannian manifolds. In our case we have to change differential forms 
and measures on Riemannian spaces into similar objects on Em 

VI. HEAT EQUATIONS FOR DISTINGUISHED TENSOR FIELDS IN VECTOR 

BUNDLES 

To generalize the results presented in Section IV to the case of d-tensor fields in £/y we 
use the Ito idea of stochastic parallel transport [27,30] (correspondingly adapted to trans- 
ports in vector bundles provided with N-connection structure). 

A. Scalarized Tensor d-fields and Heat Equations 

Consider a compact bundle E^ and the bundle of orthonormalized adapted frames on 
Ex denoted as O {Ejy) . Let {L ,L 1 , ...,L q _i} be the system of canonical horizontal vector 
fields on O (£V) (with respect to canonical d-connection T j3y . The flow of diffeomorphisms 
r {t) = r (t, r, w) on O {Em) is defined through the solution of equations 

dr(t) = L«(r (t))o6w*(t), 

r (0) = r, 

and this flow defines a diffusion process, the horizontal Brownian motion on O (En) , which 
corresponds to the differential operator 

^) = ^E4(4)' (6-i) 
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For a tensor d-field S (u) = Splfo^p* (u) we can define its scalarization F s (r) = 
(a system of smooth functions on O similarly as we have done in Section II, but in 

our case by using frames satisfying conditions (2.23) in order to deal with bundle O (£j\f) ■ 
The action of Laplace-Beltramy operator on d-tensor fields is defined as 

i. / s. \\aia 2 ...a p 

dJd.t)) =g°*t%£;:Z* 

\ / / /3l/32—/3 q 

^ 

where DT is the covariant derivation with respect to r ^. We can calculate (by putting 
formula (2.21) into (6.1) that 



A 0(^)1^5/3 1 /3 2 ..aJ - \*U3)lhlh...l} q ■ 

For a given d-tensor field S — S (u) let defined this system of functions on [0, oo) x O (£jv) : 

v$g::£ (*, r) = E [F a s iXX ir (*> < ■ 

According to the theorem 4.1 V^^ 2 "^ p is a unique solution of heat equation 

dV 1 

«■ = -^o lM V, (6.2) 

K|t=0 - *Sfcfo...p q - 

In a similar manner we can construct unique solutions of heat equations (6.2) for the 
case when instead of differential forms one considers 7£ m+n -tensors (see [7] for details con- 
cerning Riemannian manifolds). We have to take into account the torsion components of the 
canonical d-connection on £^ . 

B. Boundary Conditions 

We analyze the heat equations for differential forms on a bounded space £^ : 

da 1 . 

w = 2 Da ' (6 ' 3) 

a\t=o = f, 
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O^norm 0, {d,Ol) norm 0, (6-4) 

where □ is the de Rham-Codaira Laplacian (5.8), 

Unarm = Oq-1 (u) dll^ 1 , 

We consider the boundary of to be a manifold of dimension q = m + n and denote by 
the interior part of and as d£_\f the boundary of S^f. In the vicinity U of the boundary 
we introduce the system of local coordinates u = {(u a ) > 0} for every u G U and 

ueU n dSjsr if and only if u^ 1 = 0. 

The scalarization of 1-form a is defined as 

[F a ] p _ (r) = 6p («) 4 r = (V 3 , eg G O (£V) . 

Conditions (6.4) are satisfied if and only if 

and 

4-^=1 [*'a]f L (r) = 0, 

a=0,l,2,...,q-l, where e§ is inverse to e% 

Now we can formulate the Cauchy problem for differential 1-forms (6.3) and (6.4) as a 
corresponding problem for 7^™ +m -valued equivariant functions Vg_ (t, r) on O (£/v) : 

(t, r) = -{ Aa^) Va (t, r) + itg (r) V£ (f , r) , (6.5) 

K (0, r) = (Ff) a (r) , (P = 0, 1, g - 2), (a,£ = 0, 1, g - 1), 

(3 & (3 

where i?f (r) is the scalarization of the Ricci d-tensor and dO (£/y) = {r = (u, e) G O (£/v) , 
u G <9£V }• 
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The Cauchy problem (6.5) can be solved by using the stochastic differential equations 
for the process (U (t) ,c(t)) on 1Z n + m x U {n+m ^ : 



dV? = e%(t)o6B0(t)+SZ_ 1 6<p(t), 
del(t) = -?l(U(t))e}(t)o5U?(t) = 

(U (0) el (t) 4 (t) o SET (t) - ~T a q _ lT (U (t)) e£ (t) Sip (t) , (6.6) 

(£,f=l,2,...,<Z-l) 

where B a (t) is a (n + m)-dimensional Brownian motion, [/ (t) is a nondecreasing process 
which increase only if U {t) G dS^ . In [7] (Chapter IV, 7) it is proved that for every Borel 
probability measure /i on 7^™ +m x 7^( n + m ) 2 there is a unique solution (U (t) ,c(t)) of equations 
(6.6) with initial distribution /i. Because if 

G a/3 (C/(0)K(0)eg (0)=^ 



G^(C/(0)e- (0 eg (t) = ^a.s. 



then for every t > 

(TT (+\ \ U\ 

(this is a consequence of the metric compatibility criterions (2.10)) we obtain a diffusion 
process r (t) = (U (t) ,c(t)) on O ■ This process is called the horizontal Brownian motion 
on the bundle O (Sj^) with a reflecting bound . Let introduce the canonical horizontal fields 
(as in (2.21)) 

(v)w=^-r;( M ) ei e«, r = (11 ,e), 



T 



define the Bochner Laplacian as 

q ~ ~ 

a=l 

and put 



afl-^-i (r) = G'-^u), = -e^f q _ Xr («) G^ 1 . 
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Theorem 6.1. Let r (t) = (U (t) ,c(t)) be a horizontal Brownian motion with reflecting 
bound giving as a solution of equations (6.6). Then for every smooth function S (t,r) on 
[0,oo) x O (£j\r) we have 

dS(t,r(t))=L~S (t,r(t))6B°+ 

{3 { A o {£j ,)S) (t, r (t)) + ~^(t,r (t))}dt + (U^S) (t, r (*)) Sip (t) , 
where U q -\ is the horizontal lift of the vector field U q _i = defined as 

and 

5U q - 1 (t) 5U q ~ l (t) = a q - lq - 1 (r (*)) dt, 5U q - 1 (t) 5e% (t) = a q f ip (r (t)) dt. 

The proof of this theorem is a straightforward consequence of the Ito formula (see (a6) 
in the Appendix) and of the property that Y.aC (*) ^(t) = G a(S (U (<)) (see (2.26)). 

Finally, in this subsection, we point out that for diffusion processes we are also dealing 
with the so-called (A,L)-diffusion for bounded manifolds (see, for example, [7] and formula 
(al2)) which is defined by second order operators A and L given correspondingly on Ex and 
dEu ■ 

VII. DISCUSSION 

In the present paper we have given a geometric evidence for a generalization of stochastic 
calculus on spaces with local anisotropy . It was possible a consideration rather similar to 
that for Riemannian manifolds by using adapted to nonlinear connection lifts to tangent 
bundles [31] and restricting our analysis to the case of v-bundles provided with compatible 
N-connection, d-connection and metric structures. We emphasize that in the so-called al- 
most Hermitian model of generalized Lagrange geometry [11,12] this condition is naturally 
satisfied. As a matter of principle we can construct diffusion processes on every space E^ 
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provided with arbitrary d-connection structure. In this case we can formulate all results 
with respect to an auxiliary convenient d-connection, for instance, induced by the Christof- 
fel d-symbols (2.15), and then by using deformations of type (2.16) (or (2.17)) we shall find 
the deformed analogous of stochastic differential equations and theirs solutions. 

When the results of this paper have been communicated during the Iasi Academic Days, 
Romania, October 1994 [32] Academician R. Miron and Professor M. Anastasiei pointed our 
attention to the pioneer works on the theory of diffusion on Finsler manifolds with applica- 
tions in biology by P.L. Antonelli and T.J.Zastavniak [1,33]. Here we remark that because 
on Finsler spaces the metric in general is not compatible with connection the definition of 
stochastic processes is very sophisticated. Perhaps, the uncompatible metric and connection 
structures are more convenient for modeling of stochastic processes in biology and this is 
successfully exploited by the mentioned authors in spite of the fact that in general it is 
still unclear the possibility and manner of definition of metric relations in biology. As for 
formulation of physical models of diffusion in anisotropic media and on locally anisotropic 
spaces we have to pay a due attention to the mutual concordance of the laws of transport 
(i.e. of connections) and of metric properties of the space, which in physics plays a crucial 
role. This allows us to define the Laplace-Beltrami, gradient and divergence operators and 
in consequence to give the mathematical definition of diffusion process on la-spaces in a 
standard manner. 
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APPENDIX A: STOCHASTIC EQUATIONS AND DIFFUSION PROCESSES ON 

EUCLIDEAN SPACES 



We summarize the results necessary for considerations in the original part of this paper. 
Details on stochastic calculus and diffusion can be found, for example, in works [7,18-23]. 



Let consider the probability space (Q, J 7 , P) , where (Q, J 7 ) is a measurable space, called 
the sample space, on which a probability measure P can be placed. A stochastic process 
is a collection of random variables X = {X L ;0 < t < 00} on (fi, J 7 ) , which take values 
in a second measurable space (S, B ) , called the state space. For our purposes we suggest 
that (S, B ) is locally a r-dimensional Euclidean space equipped with a a-field of Borel sets, 
when we have the isomorphism S = W and B = B (1Z r ) , where B (U) denotes the smallest 
a-field containing all open sets of a topological space U. The index t e [0, 00) of the random 
variables X t will admit a convenient interpretation as time. 



We equip the sample space(f2, J 7 ) with a filtration, i.e. we consider a nondecreasing 
family {J 7 , t > 0} of sub a-fields of T : T s C T t Q T for 0< s < t < 00. We set J 7 ^ = 



One says that a sequence X n converges almost surely (one writes in brief a.s.) to X if 
for all weQ, excepting subsets of zero probability, one holds the convergence 

lim X n (u) =X(uj). 

n— >oo 

A random variable X t is called p-integrable if 



(X t is integrable if (al) holds for p — 1). For an integrable variable X the number 



1. Basic Concepts and Notations 





(al) 



n 
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is the mathematical expectation of X with respect to the probability measure P on (Q, JF) . 
Using a sub- cr-field Q of <r-field F we can define the value 



E(X,Q) = J X (u) du 



Q 

called as the conditional mathematical expectation of X with respect to Q. 

Smooth random processes are modeled by the set of all smooth functions w : [0, oo) 3 
t — > w (t) E 1Z r , denoted as W r = C ([0, oo) — > 1Z r ) . Set W r is complete and separable 
with respect to the metric 

p (w 1: w 2 ) = £ 2~ n [(max \ Wl (t) - w 2 (2)|) f\ 1 

n=l L ^ ' ■ 

Wi, u>2 G lU r , where a A 1 = min{a, 1}. 

Let i3 (W r ) be a topological a-field. As a Borel cylindrical set we call a set -B C W r , 
defined as B = {w : (w fa) , w (t 2 ) ,-,w (t n )) and E C B (K nr ) . We define as C C B (W r ) 
the set of all Borel cylindrical sets. 

Definition Al. A process {X t ,T t -, < t < oo} is said to be a submartingale (or 
supermartingale ) if for every < s < t < oo, we have P— a.s. E (X^J 7 ) > X s {or 
E (XtlJ 7 ) < X s ). We shall say that {X t ,jF t , < t < oo} is a martingale if it is both a 
submartingale and a supermartingale. 

Let the function p (t, x) , t > 0, x E W is defined as 

i2n 



p (t, x) = (2nt) 2 exp 



\x\ 



2t 



and X = (X t ) te [ ^ is a r-dimensional process that for all < t± < ... < t m and Ei E B (1Z r ) , 
i = 1,2,..., ra, 

P{X tl E E 1 ,X t2 E E 2 , ...,X tm E E m } = 
J (j, (dx) J pfa,xi- x) dxi J p(t 2 -t 1 ,x 2 - xi) dx 2 ... J p(t m - * m _i, x m - x m _i) dx m , 



(a2) 
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where \x is the probability measure on(7?. r , B (W)) . 

Definition A2. A process X = (X t ) with the above stated properties is called a r- 
dimensional Brownian motion (or a Wiener process with initial distribution /i. A probability 
P x on (W r ,B (W r )), where P{w : w (h) E E 1: w(t 2 ) E E 2 ,...,w(t m ) E E rn } is given by 
the right part of (a. 2) is called a r-dimensional Wiener distribution with initial distribution 
fi. 

Now, let us suggest that on the probability space (fi, J 7 , P) a filtration (jF t ) , t E [0, oo ) 
is given. We can introduce a r-dimensional (jF t )-Brownian motion as a d-dimensional smooth 
process X = (X (0)tg[ooo) > (-^-coordinated and satisfying condition 



£[exp [i<£,X t -X s >] \F S 



cxp 



(t-s) |e| 2 /2 



a.s. 



for every £ E lZ r and < s < t. 

The next step is the definition of the Ito stochastic integral [21,7,18-20]. Let denote as 
£ 2 the space of all real measurable processes $ = {$ (t, u)} t >o on fi, (!F t ) -adapted for every 
T > 0, 



$ ||;> = E 



1 

j $ 2 (s,u) ds < oo 



where "= " means "is defined to be". For $ E £2 we write 



*I| 2 =E 2 " 



n=l 



2,n 



AO- 



We restrict our considerations to processes of type 

00 

$ (t, u) = f (u) I {t=0} M (*) > 



(a3) 



where J A (fi) = 1, if A C B and J A (5) = 0, if A C £. 

Let denote .M2 = {-X" = (X t ) t>0 ;X is a quadratic integrable martingale on (fi,JF, P) 
referring to (J 7 t) t >o anc ^ -Xo = a.s. } and M. c 2 = { X E M.2', t — > X is smooth a.s.}. For 
X E M.2 we use denotations 



= E 



X 



,T>0, 
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and \X\ = £ 2-»(|X| n Al). 
n=i 

Now we can define stochastic integral on (jF t )-Brownian motion B (t) on (Q, JF, P) as a 
map 

For a process of type (a3) we postulate 

n-1 

/ ($) (t, W ) = Yl ifi M ( 5 - 5 ")) + /n M (*, ™) - B (t n , CO))] = 

oo 

E/*( s (*A*i+i)- s (*A*i)) ( a4 ) 



i=0 



i=0 



for t n <t < t n+1 ,n = 0, 1, 2, .... 

Process J ($) G defined by (a4) is called the stochastic integral of $ G £2 on 
Brownian motion B (t) and is denoted as 



I ($) (f) = J $ ( s , w ) d£ ( S) w )= J <$> (s) dB (s) . 



(a5) 



It is easy to verify that the integral (a5) satisfies properties: 



|J($) 



2,T- 



and 



00 ,. 
E(l($)(t) 2 )=Y,E[f?(tf\t l+1 -tf\t l )]=E U 2 (s,w)ds 

i=0 n 



/ (a$ + /TO) (0 = al ($) (t) + pi (*) (f) , 



for every $, * G £ 2 (a, (3 E K) and t > 0. 

Consider a measurable space JF) equipped with a filtration ) . A random time 
T is a stopping time of this filtration if the event {T < t} belongs to the <T-field(jF t ) for 
every t > 0. A random time is an optional time of the given filtration if {T < t} G (Tt ) 
for every t > 0. A real random process X = (X t ) t>0 on (0, : T, P) is called a local (^i) 
-martingale if it is adapted to (Tt ) and there is a sequence of stopping (Tt ) -moments o n 
with a n < 00, a n | 00 and X n = (X n (t)) is a )-martingale for every n — 1,2, where 
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X n — X (t A cr n ) . If X n is a quadratic integrable martingale for every n, than X is called a 
local quadratic integrable [T t )-martingale . 

Let denote A4 l 2 oc = { X = (X t ) t>Q ,X is a locally quadratic integrable [Tt )-martingale 
and Xq = a.s.jand M. C 2° C ={X G M l 2 oc : i — > X t is smooth a.s.}. In a similar manner 
with the Brownian motion we can define stochastic integrals for processes $ G C 2 and 
$ G £3°° on M C Ai l 2° (we have to introduce M (tj,cu) instead of B (tj,u) respectively for 
tj = t i+ i,tj = ti,ti = t n ,tj = 1 in formulas (a4)). In this case one denotes the stochastic 
integral as 

t 

J M ($) (f) = J $( s )dM(s). 


A great part of random processes can be expressed as a sum of a mean motion and 
fluctuations (Ito processes) 

t t 
X(t)=X(0) + J f{s) ds + Jg(s)dB(s), 


t t 
where / / (s) ds defines a mean motion , f g (s) dB (s) defines fluctuations and / dB is a 


stochastic integral on Brownian motion B (t) . In general such processes are sums of processes 
with limited variations and martingales. Here we consider the so-called smooth semimartin- 
gales introduced on a probability space with a given filtration {T t ) t>Q , M i (t) G Mj oc and 
A % {t) being smooth (jF t ) -adapted processes with trajectories having a limited variation and 
A l (0) = (i=l,2,...,r). So a smooth r-dimensional semimartingale can be written as 

X i (t) = X i (0) + M l (t) + A 1 (t) . (a6) 

For processes of type (a6) one holds the Ito formula [21,7,18-20] which gives us a differential- 
integral calculus for paths of random processes: 

t t 
F (X (*)) — F (X (0)) =Jf;(X (s)) dM l (s) + - J F» (X (s)) d < MAI' > (s, ) , (a7) 



where F( = §^,F"j = d ^Q x j , < M*M J > is the quadratic covariation of processes M % 

,Mi G M.2-, which really represents a random process A = A t , parametrized as a difference 
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of two natural integrable processes with the property that M t N t — A t is a (.F^-martingale. 
Here we remark that a process Q = Q t is an increasing integrable process if it is {T t )- 
adapted, Q = 0, the map t — > Q\ is left-continuous , Q t > and £ (Qt) < oo for every 
i G [0,oo). A process Q is called natural if for every bounded martingale M = (M t ) and 
every t G [0, 00) 





T 




- t 


E 


J M s dA s 


= E 


J M s dA s 




.0 




.0 



There is another way of definition of stochastic integration instead of Ito integral , the 
so-called Fisk-Stratonovich integral, which obeys the usual rules of mathematical analysis. 
Let introduce denotations : A is the set of all such smooth (jF t ) -adapted processes A = (A t ) 
that Aq = and t — > A t is a function with limited variation on every finite integral a.s.; 
B is the set of all such (jF t )-predictable processes $ = (<& t ) that with the probability one 
the function t — > $ t is a bounded function on every finite interval and (t, a;) — > X t (a;) is 
C/-B (7?. r )-measurable; (9 is the set of quasimartingales (for every X C (9 we have the 
martingale part Mx and the part with limited variation ). For every $ G B and X G (9 one 
defines the scalar product 

t t 
($ o X) = X (0) + J $ (s, iu) rfM x (s) + y $ (s, (s) , t > 0, 



as an element of O . One introduces an element $ o dX G dO as 

$dX = $ o dX = d ($ o X) 
in order to define the symmetric Q-product : 

Y o dX = YdX + \dXdY 

t 

for dX G dO and F 6 0. The stochastic integral JY o dX is called the symmetric Fisk- 



Stratonovich integral. 
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2. Stochastic Differential Equations 



Let denote as A c ' r the set of functions satisfying the conditions : a (t, w) : [0, oo) x W r — > 
TZ r x TZ r are B ([0, oo)) x B {W r ) /B {TZ r <g> ft r )-measurable, for every t G [0, oo) a function 
W r 3 w -> a (t, w) G TZ r x ft r is B t (W r ) /B {ll d <g> ft c ) -measurable, where W x TZ C denotes 
the set of r x c matrices, B (1Z r <g> 1Z r ) is the topological a-field on W x 1Z C , obtained in the 
result of identification of 1Z r x 1Z C with rc-dimensional Euclidean space. 

Suppose that values a G A r ' c and /3 G A r,x are given and consider the next stochastic 
differential equations for a r-dimensional smooth process X = (X (t)) t>Q : 

r 

dX e t =J2 a i (*> x ) rf57 (*) + £ (*> x ) dt ' ( a8 ) 
j'=i 

(e= 1,2, ...,r), 

for simplicity also written as 

dX t = a (t, X) dB (t) + (3 (t, X) dt. 

As a weak solution (with respect to a c-dimensional Brownian motion B (t) ,B (0) = 
a.s.) of the equations (a8) we mean a r-dimensional smooth process X = (X (t)) t>0 , defined 
on the probability space T ,P) with such a filtration of cr-algebras (.Fu) u>/ that X — X (t) 
is adapted to (J r u ) u>/ , i.e. a map G f2 — > X(cj) G W r is defined and for every t G [0, oo) 
this map is TujB u (W r )-measurable; we can define processes <3>jg (t, a;) = (t, X (u)) C £3°° 
and # 5 (t, w) = (3 s (t, X (u)) C £i° c ; values X (i) = (X 1 (t) , X 2 (t) , X r (t)) and B (t) = 
(B 1 (t) , B 2 (t) , B c (*)) satisfy equations 

X* (t) - X* (0) = ^ / (s, X) dB" (s) + f f3 s (s, X) ds, (a9) 

o 

with the unit probability, where the integral on dB 13 (s) is considered as the Ito integral (a7). 
The first and the second terms in (a9) are called correspondingly as the martingale and drift 
terms. 

Let a (t, x) = <jj (t, x) be a Borel function (t, x) G [0, oo) x lZ r — > lZ r <E> TV and b (t, x) = 
{V (t, x)} be a Borel function (t, x) G [0, oo) x ft r -> W . Then a (t, w) = a (t, w (t)) C „4 r,c 
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and (3 (t, w) = b (t, w (t)) G A 1 "' 1 . In this case the stochastic differential equations (a8) are of 
the Markov type and can be written in the form 

r 

dx i (t) = Y, a i (*> x (*)) dBk (*) + &i (*> x (*)) ( al °) 

fc=l 

If a and 6 depend only on x G we obtain a equation with homogeneous in time t 
coefficients. 

Function $ (x, w) : 7Z r xW c -> W^ r , W C = G C ([0, oo) -> ft r ) ; w (0) = 0} is called 
£ (7£ r x W / c )-measurable if for every Borel probability measure fx on W there is a func- 
tion ^(x,™) : ft r xW£ -> W c , which is x W C ) M /B (W r )- m easurable for all 

x (/i) , $ (x, u>) = $ M (x, u>) and P^-almost all w (P w is the c-dimensional Wiener measure 
on Wq, i.e. a distribution 5 ). 

A solution X = (X (£)) of the equations (a8) with a Brownian motion B — B (t) is called 
a strong solution if there is a function F(x,w) : VJ xW c -> W r , which is £(ft v x W C )- 

pW 

measurable for every x & 7Z r , w —> F (x,w) is B t (Wq) /B t (W c )-measurable for every 
t > and X = F (X (0) , B) a.s. 

We obtain a unique strong solution if for every r-dimensional(jF u )-Brownian motion 
B = B (t) (B (0) = 0) on the probability space with the filtration (F t ) and arbitrary (Fq)- 
measurable 7£ r -valued random vector X — F (£, B) is a solution of (a8) on this space with 
X (0) = £ a.s. So, a strong solution can be considered as a function F (x, w) which generates 
a solution X of equation (a8) if and only if we shall fix the initial value X (0) and Brownian 
motion B. 

3. Diffusion Processes 

As the diffusion processes one names the class of processes which are characterized by 
the Markov property and smooth paths (see details in [22, 23, .7]). Here we restrict ourselves 
with the definition of diffusion processes generated by second order differential operators on 
U r : 
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Af (x) = \ ± a- (x) ^LL. ( x ) + ± U (x) 9 J- (x) , (all) 

where a lJ (x) and b l (x) are real smooth functions on TV, matrix a lJ (x) is symmetric and 
positively defined . Let denote by TV = TV U{A} the point compactification of TV. Every 
function / on TV is considered as a function on TV with / (A) = 0. The region of definition 
of the operator (all) is taken the set of doubly differentiable functions with compact carrier, 
denoted as C\ (TV) . Let B (w r ) be the a-field generated by the Borel cylindrical sets, where 
W r = {w : [0, oo) 3 t -> w (t) G TV is smooth and if w (t) = A , then w (?) = A for all 
t' > t. The value e (w) = inf{t; w (t) — A, w e W r } is called the explosion time of the path 
w. 

Definition A3. A system of Markov probability distributions {P x ,x e TV} on 
(W r ,B(W r ^, which satisfy conditions : Px{w : w(0) — x} — 1 for every x G VJ'\ 
f (w (t))-f(w (0)) -J (Af) (w (s)) ds is a (p x , B t (W r ))-martmgale for every f e C\ (TZ r ) 
and x G lZ r defines a diffusion measure generated by an operator A (or A- diffusion) . 

Definition A4. A random process X = (X (t)) on TV is said to be a diffusion process, 
generated by the operator A (or simply a A- diffusion process) if almost all paths [t — > X (t)] 
G W r and probability law of the process X coincides with (■) — f P x (■) jj, (dx) , where 
H is the diffusion measure generated by the operator A and {P x } is the probability law of 
X(0). 

To a given A-diffusion we can associate a corresponding stochastic differential equation. 

Let the matrix function a (x) = (a) (x)) G TV x TV defines a lj (x) — ( x ) a { ( x ) anc ^ 

v J ' k=i 

consider the equations 

r 

dx i (t) = <?i ( x (*)) dBk (*) + u ( x (*)) dt ( a12 ) 

k=l 

There is an extension of (Q, T,P) with a filtration (F t ) of the probability space 

(\V r ,B(w r ),P x ) 

and with a filtration B t (W r ^j and a (jF t )-Brownian motion B (t) (see [7] and the previous 
subsections in this Appendix) that putting X (t) = w (t) and e = e (w) one obtains for 
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t E [0,e) 

r t t 
X 1 (t) = x t + J2 J 4 (X (s)) dB k (s) +fv(X (s)) ds. 

k=1 

So (X (t) ,B(t)) is the solution of the equations (al2) with X (0) = x. 

If bounded regions are considered, diffusion is described by second order partial differen- 
tial operators with boundary conditions. Let denote D = VJ t = {x = ( 0C j X j .... Ob r )-x r > 
}, 3D = {x e D,x r = 0},D° = {x e D;x r > 0}. The Wentzell bound operator is defined 
as a map from C\ (L) to the space of smooth functions on dD of this type: 

Lf (x) = \ £ « y (x) (x) + £ (*) ^ (x) + v (x) ^ (x) - p (x) Af (x) , (a!3) 

i , j — 1 2 — 1 

where x G <9.D, a*- 7 (x) , /3 l (x) , p (x) and p (x) are bounded smooth functions on dD, a %: > (x) 
is a symmetric and nondegenerate matrix, p (x) > and p (x) > 0. 

A diffusion process defined by the operators (all) and (al3) is called a (A, L) -diffusion. 
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